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Abstract 

In the present paper we show how non-classical, quite accurate, critical 

exponents can be extracted in a very simple way from the Pade analysis 

of the results obtained by mean field like approximation schemes, and in 

particular by the cluster variation method. We study the critical behavior of 

the Ising model on several lattices (quadratic, triangular, simple cubic and 

face centered cubic) and two problems of surface critical behaviour. Both 

unbiased and biased approximants are used, and results are in very good 

agreement with the exact or numerical ones. 
PACS numbers: 05.50.+q 
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As is well known, mean field like approximations are very useful tools in the inves- 
tigation of phase transitions, but they fail completely in predicting critical exponents of 
low-dimensional systems, giving results which are independent of the model and the dimen- 
sionality. This is due to neglecting long range correlations, and thus they are completely 
unreliable if one is interested in the behavior of a model near critical or multicritical points. 

On the other hand, there are some mean field like approximations which describe very 
accurately the low- and/or high-temperature behaviour of statistical mechanical models. In 
particular, the Cluster Variation Method (CVM) introduced by Kikuchi yj and reformulated 
several times has been shown || to reproduce exactly many terms of the high- and low- 
temperature expansions of thermodynamical quantities like specific heat, magnetization and 
susceptibility. 

In the present paper we show how this property of the CVM (and also of other schemes, 
as we shall see) can be used to build up a very simple procedure for determining quite 
accurately the critical temperature and critical exponents of a given model, which relies on 
the Pade analysis of low- and high-temperature results obtained by the CVM. The plan 
of our presentation is as follows: we first review the main ideas of the CVM and of Pade 
approximants, then we describe in detail our technique and several test applications and 
finally we discuss results and possible generalizations. 

The Cluster Variation Method in its modern formulation |J can be seen as a truncation 
of a cumulant expansion of the functional to be minimized in the variational principle of 
statistical mechanics. The latter states that the free energy F of a model system described 
by the hamiltonian Hona lattice A can be obtained by minimizing the functional 

F[p A ] = Tt(paH + k B Tp A \np A ), (1) 

where k B and T are as customary Boltzmann's constant and absolute temperature, with 
respect to the density matrix p A , subject to the constraint TrpA = 1- Upon introducing 
the cluster density matrices p a = Tr A \ a p A , where a is a cluster of n a sites and the trace 
is performed over all variables out of a, Eq. ([!]) is approximated by a restricted variational 
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principle for the functional 

F[{p a ,a e P}] = Tr(p a /i a ) + k B T a a Tr(p a \n p a ), (2) 

where P is a set of "maximal" clusters and all their subclusters, h a is the n a -body interaction 
contribution due to the cluster a (maximal clusters should be taken large enough to contain 
all kind of interactions present in Ti), the coefficients a a obey |§ 

a a = l, V/?eP (3) 

fiCatP 

and the p a 's are subject to the constraints 

Trp a = 1,V«GP and p a = Trp^pp, Va C (3 E P. (4) 

Local minima of F can be obtained in a simple way by means of an iteration scheme devised 
by Kikuchi 0], the Natural Iteration Method (NIM) which, for maximal clusters up to 8-9 
sites as those we used, does not require large amounts of CPU times. 

It should be clear that an approximation in this scheme is uniquely defined by the choice 
of the maximal clusters (so if the maximal clusters are all the elementary cells of a simple 
cubic lattice, we will speak of the cube approximation, and so on), and it will not be accurate 
near critical points, where the correlation length of the system becomes larger than the size 
of the maximal clusters. Generally speaking, taking larger maximal clusters will narrow the 
region in which the approximation is not accurate. 

Let us now turn to a brief review of some basic facts about Pade approximants 0. An 
[L, M] Pade approximant to a function F(z) is defined by 

\t M] = = Po+Pi z + P2 z2 + • • • +PlZ L . , 

L ' l ~Q M (z) l + qiZ + q2Z 2 + ... + qMZ M- U 

The usefulness of these approximants in statistical mechanics derives from the fact that a 
function with a power-law singular behaviour 

F(z)*A(z)(l--f) \ z^z;, (6) 



with A(z) analytic at z c , will have a logarithmic derivative of the form 

D(z) = ^-hxF{z) « — ^-[1 + 0(z - z c )}, z -> z~ (7) 

and a simple pole like that in Eq. (^) can be represented exactly by Pade approximants. 
Clearly, z c and —A are given respectively by the pole and the corresponding residue of the 
Pade approximant. Furthermore, if the exact value (and this is sometimes the case) or a very 
accurate estimate is available for z c , better estimates of A can be obtained by approximating 
the function 

\*(z) = (z c - z)D(z) n\ + 0(z-z c ), z^ z-. (8) 

In this case A is estimated directly as the value of the approximant at z c . If no estimate 
of z c is given as an input one speaks of unbiased approximants, otherwise the approximants 
are called biased. 

We can now turn to the description of our technique. To fix ideas, let us consider the 
Ising model on a face centered cubic (fee) lattice. The hamiltonian is 
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k B T 



-J^SiSj -h^Si, (9) 
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where J and h are the (reduced) interaction strength and magnetic field, Sj = ±1 is the 
z-component of a spin 1/2 operator at the lattice site i and the first summation is over 
nearest neighbours (n.n.). Our first step consists in choosing two CVM approximations 
(i.e. two sets of maximal clusters) Mi and M 2 in such a way that M 2 can be regarded 
as an improvement with respect to Mi. To be concrete, in the following M\ will be the 
octahedron plus tetrahedron approximation proposed by Aggarwal and Tanaka 0, while 
Mi will be the newly developed oriented rombohedron approximation (to be described in 
detail in a separate paper |J). This approximation is obtained by selecting as maximal 
clusters all the primitive rombohedral unit cells of an fee lattice with a given orientation 
(an fee lattice can be decomposed into primitive cells of this shape in four different ways, 
corresponding to different orientations). Since a rombohedron is made up of an octahedron 



with two tetrahedra attached on opposite faces, we can expect this new approximation to 
improve with respect to M±. Indeed, this can be easily verified by comparing different 
estimates for the critical temperature: one obtains T c = 1/J C — 10.03 from the CVM 
tetrahedron approximation, 10.01 from Mi, 9.97 from M 2 and 9.83 from high-temperature 
expansions ]7j|. 

The next step is to compare results from Mi and M2 to determine the temperature 
(or, equivalently, interaction strength) range in which Mi is accurate in some sense. Let 
us consider first the low-temperature region, in order to obtain estimates for the critical 
temperature and the critical exponent (3 associated to the vanishing of the order parameter 
m = (sj). We have chosen as a measure of the accuracy of Mi the quantity 5m(J) = 
I mi (J) — 777.2 (J) I, where m&(J) is the value of the order parameter (in zero field) as given by 
approximation Mjfc, which, not too close to the critical point, should be a good approximation 
to the absolute error of mi (J) with respect to the exact (unknown) value (it can be checked 
on two dimensional problems that this is a very reasonable assumption). Then we define a 
value J min of the interaction strength such that 8m(J) < e for J > J mm , where e is a small 
positive number, and say that Mi (and consequently also M 2 ) is accurate for J > J m i n . Of 
course it would be desirable to take e very small, since large values of e cause poor values of 
the order parameter to be treated as accurate, but, on the other hand, reducing e narrows 
the temperature range on which the Pade analysis will be made. A good compromise, which 
we have used throughout this paper is e = 10~ 5 , which for M 1 and M 2 as above yields 

Jmin 0.14. 

We then determine Pade approximants for the logarithmic derivative of the magnetization 

D(z) = j-lnm(«), (10) 

where (as customary for low temperature approaches) z = e~ J is a variable which vanishes 
at zero temperature and m(z) is given by approximation Mi (or M 2 , if this is not very 
time-consuming), by requiring that, for a given pair of positive integers L and M, 

[L,M](z n ) = D(z n ), z n = exp(J mhl + nAJ), n = 0, 1, 2, . . . L + M. (11) 
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The choice of A J is constrained, since the above linear system becomes ill-conditioned as 
the interpolation points become too close, either in the z direction (A J small) or in the D(z) 
direction (A J large). In Tab. | we report estimates of T c and (3 obtained by unbiased [L, M] 
approximants for AJ = 0.03. They indicate clearly that T c ~ 9.79 (a strong improvement 
with respect to the M 2 estimate) and (3 ~ 0.31, which are in good agreement with the high- 
temperature expansion critical temperature and with (3 = 0.3258 ± 0.0044 (Monte Carlo 
simulations for the simple cubic lattice, ref. ||) respectively. Different values of A J in the 
neighborhood of 0.03 give similar results. 

Once an estimate for T c has been obtained, it can be used to construct biased approxi- 
mants for the logarithmic derivative of the high temperature susceptibility. In the disordered 
phase, we measure the accuracy of CVM approximations by means of <5c(J) = \c%(J) — c 2 (J)|, 
where c&(J) is the n.n. correlation function (s,-Sj) as given by approximation Requiring 
<5c(J) < 10~ 5 for J < J max we obtained J max = 0.048. [L, M] approximants have then been 
determined for the function 

7 *(w) = (w c - w )-^x{w), (12) 

where w = tanh(J) and xi w ) is the uniform zero field magnetic susceptibility. The L + M + l 
interpolation points have been defined by w n = tanh( J max — nAJ),n = 0, 1, . . . L + M. For 
A J = 0.002, all the [L, M] approximants with 2<L<5andL-l<M< M+l, except the 
[2, 1] one, give 7 = 1.26, a result which is not far from the best estimate 7 = 1.2390 ±0.0071 
(Monte Carlo simulations for the simple cubic lattice, ref. 0). A similar analysis on the 
specific heat gives some indication for a ~ 0.14, but data do not accumulate very well. 
It should also be noted that biasing approximants with the high temperature expansion 
estimate for the critical temperature improves results for the critical exponents: one obtains 
(3 ~ 0.33 and 7 ~ 1.24. 

We have done a similar analysis for the simple cubic lattice, choosing for M\ the well- 
known cube approximation and for M2 a newly developed (again to be described in a separate 
paper |J) approximation which uses as maximal clusters both the elementary cubes and the 
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"stars" formed by one site surrounded by its six nearest neighbors, which is a straightforward 
generalization of the approximation used by Finel and de Fontaine |9j in their investigation 
of the two dimensional ANNNI model. J m i n and J max are 0.28 and 0.13 respectively. Low 
temperature unbiased approximants indicate clearly T c ~ 4.51, in very good agreement with 
T c = 4.511424 ± 0.000053 (Monte Carlo simulations, ref. §), while (3 is between 0.30 and 
0.31. Biased approximants for /3, however, are in favor of (3 ~ 0.31. Results from biased 
high temperature approximants for the susceptibility suggest 7 ~ 1.24. 

In Tab. |T| we report results for two dimensional lattices. In this case, due to the reduced 
dimensionality, the CVM is expected to be less accurate, but biasing the approximants 
with the exactly known critical temperature yields indeed very good results. We used 



the approximation B 2 n proposed by Kikuchi and Brush |IJ for the square lattice and a 
straightforward generalization of it for the triangular lattice, choosing N = 3 for Mi and 
iV = 4 for M 2 . 

Finally, we want to discuss two test applications of our method to surface problems. 
We have considered a semi-infinite Ising model with a (100) free surface and unmodified 



surface coupling, which is known |TT| to exhibit a so-called ordinary transition, with the 



surface layer magnetization mi vanishing with an exponent (3i which differs from the bulk 



exponent (3 and is estimated to be 0.78 ± 0.02 by Monte Carlo simulations |L2| and 0.816 



by second order e-expansion |fL3|| . For this problem we have developed another new CVM 
approximation (again described in detail in ref. |J), which we refer to as 4 x N, in which 
the semi-infinite system is approximated by a film of N layers, with the topmost layer 
representing the free surface, and the bottom layer constrained to the bulk, which in turn 
is studied in the cube approximation. In this system the maximal clusters for the CVM are 
chosen as those clusters with 4N sites, formed by a column of N — 1 elementary cubes. We 
used N = 4 for Mi and N = 5 for M 2 , obtaining J min = 0.30. Approximants biased with 
the previously determined bulk critical temperature for the simple cubic lattice indicate 
f3i ~ 0.78, in perfect agreement with Monte Carlo simulations. 

Our last test application aims to illustrate that the CVM is not the unique classical 



approximation on which our method can be based. In particular we have considered an 
approximation recently proposed by Lipowski and Suzuki |nj for two dimensional systems 
(referred to as the LS approximation), which has been shown in ref. ]15[ to yield the exact 
boundary magnetization of the square lattice Ising model. It is essentially a transfer matrix 
mean field approximation, where the boundary magnetization is calculated by considering 
two strips of width N and N — 1 with periodic boundary conditions along the infinite 
direction and applying the same effective field h c s at one side of each strip, where h e s is 
determined in such a way that the boundary magnetizations on the opposite sides of the 
strips be equal. The approximation is no more exact for more complicated models, but one 
can expect it to be quite accurate. We have applied it to the three-state Potts model in two 
dimensions, described by the hamiltonian 

n 



k B T 



-*E<W ( 13 ) 

(ij) 



We have taken N = 4 for Mi and N = 5 for Mi {K m \ n = 0.90), and have estimated 
the boundary magnetization exponent (3\ using approximants biased with the exact critical 
temperature. Our results indicate clearly (3\ ~ 0.55, in very good agreement with f3\ = 5/9, 
obtained combining Cardy's result f3i = u/(3i/ — 1) JTB| with the conjecture v = 5/6, 
supported by many numerical results Jl] 



We have shown how quite accurate estimates of critical temperatures and critical expo- 
nents can be obtained by a Pade analysis of high- or low-temperature results of mean field 
like approximations, and especially of the Cluster Variation Method. The method is quite 
simple (some analytical work can be needed to construct new CVM approximations, when 
necessary) and not at all time-consuming (all the calculations reported in this paper took 
a few hours of CPU time on a DEC Alpha machine), but nevertheless it yielded very satis- 
factory results in several test applications, also in surface problems. There are also several 
possible future developments, among which the calculation of critical amplitudes and the 
use of more sophisticated approximants like differential and partial differential approximants 
(the latter applying to the study of multicritical phenomena) are worth mentioning. 
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TABLES 



TABLE I. (T c ,/3) for the fee Ising model (unbiased approximants, J, 


nin = 0.14, 

inn " 


AJ = 0.03). 


L [L,L-1] 


[L,L] 




[Li, Li + 1J 


4 (9.75, 0.30) 


(9.79, 0.31) 




(Q 7S 30) 


5 (9.78, 0.31) 


(9.78, 0.31) 




(9.79, 0.31) 


6 (9.78, 0.31) 


(9.79, 0.31) 




(9.79, 0.31) 


7 (9.78, 0.31) 


(9.79, 0.31) 




(9.79, 0.31) 


8 (9.79, 0.31) 


(9.79, 0.31) 




(9.79, 0.31) 


TABLE II. Results for two dimensional lattices (biased approximants, exact T c ). 


Lattice «/ m in 


P 


Jm&x 


7 


Square 0.53 


0.123 


0.30 


1.73 


Triang. 0.35 


0.125 


0.17 


1.74 
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